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A SIMPLIFIED  METHODOLOGY  FOR  COMPUTING 
CHEMICAL  HAZARD  DISTANCES 


I 


1.  SCOPE 

1 . 1 Objective.  It  is  the  objective  of  this  task  to  develop  a methodology  for  the  solution 
of  the  Gaussian  diffusion  model  which  is  simple  enough  to  be  performed  on  a pocket  calculator.  It 
is  intended  that  this  would  include  the  effect  of  the  bounded  mixing  layer’  and  the  2-minute 
dosage  correction^  as  employed  with  the  G and  V agents. 

1.2  Background.  The  Methodology  for  Chemical  Hazard  Predictions,  which  was  agreed 
upon  by  the  Technical  Advisor  to  the  Chemical  Standards  Working  Group  of  the  Department  of 
Defense  Explosives  Safety  Board  in  1974,  was  published  as  DODESB  Technical  Paper  No.  10^  in 
March  1975.  The  mathematical  representation  of  this  methodology  was  stated  in  very  general  form 
and  its  inherent  complexity  prevented  its  immediate  adoption  for  field  work. 

In  1976,  a system  of  graphic  aids  was  prepared  by  the  Systems  Assessment  Office, 
Chemical  Systems  Laboratory,  and  was  published  by  DARCOM  in  March,  1977,  as  a Handbook  for 
Chemical  Hazard  Prediction.’  In  June  1977,  ARCSL-TR-77049'*  was  published  by  Systems 
Assessment  Office  to  document  a complete  computer  program  of  this  methodology. 

Since  the  Handbook’  was  developed  as  the  basic  implementing  document  of  this 
methodology,  it  was  presented  in  generalized  form  which  attempted  to  portray  all  the  variables  of 
the  model  over  the  ranges  of  interest.  This  document  was  large,  containing  some  170  graphs,  but 
still  treated  many  variables  as  factor  corrections  to  a first  estimate. 

The  computer  program,  documented  as  ARCSL-TR-77049,  was  also  developed  as  a 
basic  tool  and  attempted  to  cover  all  of  the  methodology  in  the  most  concise  manner.  As  a result, 
the  program  is  made  up  of  eight  subroutines  containing  some  1000  Fortran  statements.  The 
program  is  fast,  relatively  easy  to  use,  and  will  be  maintained  on  the  ARRADCOM  Univac  1 1 08 
computer  at  Aberdeen  Proving  Ground,  Maryland,  for  remote  use  through  the  telephone. 

There  remained,  however,  an  area  of  quick  field  applications  which  neither  of  these 
approaches  satisfied  completely.  An  approach  was  sought  which  would  approximate  some  of  the 
more  complex  procedures  in  the  methodology  and  lead  to  simple  statements  which  could  be 
programmed  on  the  pocket  calculator. 

1.3  Approach.  Once  the  complete  methodology  was  available  on  a high-speed 
computer,  investigation  began  as  to  how  some  of  the  more  complex  operations  could  be 
approximated.  The  transition  from  the  Gaussian  to  the  Box  Model*  was  of  prime  concern  since  this 
had  not  been  fully  covered  by  the  graphic  solutions  presented  in  the  handbook.  From  the  many 
curves  in  the  handbook,  it  was  evident  that  this  transition  occurred  over  a relatively  short  distance 
and  an  approach  was  sought  to  define  this  zone  as  a function  of  stability  and  height  of  the  mixing 
layer.  Once  the  transition  zone  was  defined,  the  dosage-versus-distance  curve  would  be 
approximated  by  three  line  segments  representing  the  Gaussian  Model,  mixing-layer  transition,  and 
Box  Model. 


Diffusion  under  an  inversion  cap  with  vertical  distribution  uniform. 
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In  order  to  simplify  the  2-minute-coirected  methodology  as  developed  for  ORG  40^ 
and  reproduced  in  Technical  Paper  No.  1 0,^  the  full  iterative  program  was  executed  and  the  final 
dosage  correction  factor  at  each  distance  was  converted  to  an  effective  time  of  exposure.  This 
effective  time  was  then  fitted  as  a function  of  stability,  windspeed,  and  distance,  and  a simple 
algebraic  statement  of  the  2-minute  corrections  was  substituted  into  the  diffusion  equation  and 
solved  directly  for  distance. 

Finally,  the  mixing  layer  and  2-minute  corrections  were  combined  to  complete  the 
model.  Simple  algebraic  expressions  were  used  to  approximate  these  unknown  functions  with  the 
condition  that  deviations  would  be  biased  to  produce  a safe-sided  estimate  of  the  hazard  distance. 

2.  METHODOLOGY 


2. 1 Effect  of  the  Mixing  Layer.  When  a Gaussian  cloud  encounters  a physical  boundary 
such  as  the  ground  surface  or  an  elevated-inversion  cap,*  the  fraction  of  the  cloud  reaching  the 
boundary  is  reflected  back  to  reinforce  the  concentration  and  increase  the  dosage  within  this  mixing 
layer.  A model  based  on  multiple  reflections  from  these  bounding  surfaces  was  outlined  in 
ORG  17^  in  1958.  This  approach  was  reproduced  in  Technical  Paper  No.  10^  indicating  that  the 
factor  which  augments  the  dosage  at  any  downwind  distance  can  be  expressed  as  follows: 


1 

Fr  = — 2 

R 2i=i 


exp 


1 /2iH,„-H-Z 

2 \ <*z 


+ exp 


J_ 

2 


/2iHm-H-f  Z 


+ exp 


r , .2-1 

1 \ 

+ exp 

1 /2iHn,  + H + z\ 

- ^ \ ‘'Z  / _ 

. H / J 

(1) 


where 


Fj^  is  the  reflective  contribution  to  the  dosage  at  distance,  x 
is  the  height  of  the  mixing  layer  (m) 

H is  the  height  of  the  source  (m) 

Z is  the  height  of  the  sampler  (m) 

/ ^ P 


= «*zr 


, the  vertical  distribution  of  the  cloud  as  defined  in  the  handbook*  (m) 


When  one  is  interested  in  dosages  near  the  ground,  Z can  be  taken  as  zero  and  the 
expression  simplifies  to  the  following: 


Fr  = 2 

^ i=l 


/ 2iH^f  H 
\ 


2 


(2) 
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With  reasonable  values  of  H and  this  expression  equals  zero  for  small  values  of  x 
and  then  increases  with  x,  approaching: 


From  this  information,  the  beginning  of  the  transition  zone,  x j , can  be  defined  as 
the  value  of  x corresponding  to  Fr  = 0 + AE,  where  AE  is  some  specified  small  difference.  In  like 
manner,  the  extent  of  the  transition  zone,  X2,  can  be  defined  as  the  value  of  x corresponding  to 


Fr=  /£ 


2 H, 


- 1 +AE.  This  limit  is  also  the  distance  where  the  Gaussian  Model,  with  reflection. 


m 


differs  from  the  Box  Model  by  AE. 


A computer  program  was  written  to  select  values  of  Xj  and  different 

stabilities  and  heights  of  the  mixing  layer.  These  results  are  shown  in  tables  1 and  2 for  AE  = 0.01. 


TABLE  1.  VALUE  OF  xj  AS  A FUNCTION  OF  STABILITY  AND 


»m 

(m) 

A 

B 

C 

D 

E 

F 

100 

282 

550 

1023 

2138 

3548 

7079 

200 

468 

1096 

2188 

4786 

8511 

17783 

300 

617 

1660 

3467 

7762 

14125 

30903 

400 

759 

2188 

4786 

10965 

19953 

44668 

600 

1023 

3311 

7413 

17783 

33113 

77624 

800 

1259 

4467 

10233 

24547 

47863 

112201 

1000 

1479 

5495 

13188 

32359 

63095 

151355 

2000 

2399 

10965 

28840 

72443 

151355 

380188 
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TABLE  2.  VALUE  OF  X2  AS  A FUNCTION  OF  STABILITY  AND 


Hm 

(m) 

A 

B 

C 

D 

E 

F 

100 

427 

955 

1862 

4074 

6918 

14454 

200 

692 

1905 

4074 

9120 

16596 

36308 

300 

933 

2884 

6310 

14791 

27542 

63095 

400 

1148 

3802 

8710 

20893 

39811 

91201 

600 

1514 

5754 

13804 

33113 

66069 

158489 

800 

1862 

7586 

19055 

46773 

93325 

234422 

1000 

2188 

9550 

23988 

60256 

123026 

316226 

2000 

3631 

19055 

52481 

138038 

295119 

794324 

It  was  found  that  the  data  shown  in  tables  1 and  2 could  be  fitted  by  expressions  of 
the  following  form: 

1//3 


“ ^1 


^2  = ^2 


UP 


(3) 

(4) 


where 


Xj  and  X2  defines  the  transition  zone  (m) 
Cj  and  C2  are  constants  for  each  stability 
P is  the  slope  of  the  sigma  Z curve 


Given  equations  3 and  4,  one  may  now  define  three  segments  to  represent  the 


dosage-distance  relationship: 

Segment  1 

X < Xj 

0 

II 

Segment  2 

Xj  < X < X2 

Fj^  is  variable 

Segment  3 

X 

V 

X 

II 

1 
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The  derived  values  of  Cj  and  C2  are  listed  in  table  3 as  a function  of  stability.  The 
diffusion  parameters,  as  taken  or  derived  from  Technical  Paper  No.  10,^  are  also  tabulated  for  use 
in  the  equations  5 and  6.  The  reference  a values  have  been  transformed  to  one  meter  as  will  be 
discussed  under  equation  6.  Values  of  Oy  are  given  for  both  instantaneous  (ins)  and  continuous 
(con)  sources  as  discussed  in  the  handbooK.^ 


TABLE  3.  METEOROLOGICAL  PARAMETERS 


Parameter 

A 

B 

C 

D 

E 

F 

Oyj  (ins)* 

0.09 

0.0633 

0.048 

0.0634 

0.0754 

0.0796 

Oyi (con)** 

0.27 

0.1899 

0.125 

0.1268 

0.1508 

0.1592 

°z\ 

0.0222 

0.11 

0.119 

0.0898 

0.0879 

0.0791 

a 

1. 

1. 

1. 

0.9 

0.8 

0.7 

P 

1.4 

1. 

0.9 

0.85 

0.8 

0.75 

c, 

10.5 

5.5 

6.13 

9.49 

11.2 

15.3 

C2 

15.9 

9.55 

11.2 

18.1 

21.9 

31.1 

*ins  = instantaneous 
**con  = continuous 


2.2  Total  Dosage:  Segment  1,  x < Xj 

If  one  considers  total  dosage  near  the  ground  from  a point  source  release  near  the 
ground,  the  dosage-distance  relationship  is  defined  as  follows: 


where 


D(x)  = 


Q 

607T  Oy(x)  a^(x)  u 


D(x)  is  the  total  dosage  at  x (mg-min/m^) 
Q is  the  source  strength  (mg) 
u is  the  windspeed  (m/sec) 

Oy(x)  = Oyj.  ^ “ (m) 
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(5) 


(m) 


By  transferring  the  reference  distances,  Xyj.  and  x^j.,  to  one  meter,  this  is. further- 

simplified  to: 

D(x)  = ^ (6) 

607T  Oyi  u x“  ^ 

where 


Oyl  = ay,(Xy,r« 

^zl  ~ ^zr  ^^zr^  ^ 


or  the  inverse  solution : 


X = 


607T  Oyi  a^i  D u 


l/a+|3 


X < X| 


(7) 


where  D is  the  dosage  of  interest  (mg-min/m^). 

2.3  Total  Dosage:  Segment  3,  x > X2 

The  total  dosage  in  segment  3 can  be  calculated  from  the  following: 


or 


D(x)  = 


60./27T  Oy,  x« 


X = 


60.  /27t  Oyi  D 


1/a 


X >x- 


(8) 


(9) 


2.4  Total  Dosage:  Segment  2,  Xj  < x < X2 

Considering  that  the  distance  between  Xj  and  X2  is  relatively  small,  it  was  decided  to 
approximate  this  segment  as  a log-log  straight  line  between  points  (Dj , Xj)  and  (D2,  X2).  The  value, 
Dj,  is  taken  at  Xj  on  segment  1 and  D2  is  taken  at  X2  on  segment  3.  A straight  line  segment  will 
slightly  overestimate  the  dosage  in  this  region  and  thus  provide  a conservative  (safe-sided)  estimate 
of  the  distance. 


Employing  equations  6 and  8, 


Dl 


Q 

188.5  Oyi  u Xj“^^ 


(10) 
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'1  mi  1 1 |i  t 


then 


and 


Q 

150.4  Oyj  Hj^  u X2® 


In  (X2/X1) 

S = 

In  (D2/DJ ) 


(11) 


(12) 


(13) 


3.  TWO-MINUTE  CORRECTION 

3.1  Instantaneous  Sources  with  2-Minute  Correction.  A copy  of  the  program  for 
computing  dosage  with  the  2-minute  correction  for  G and  V agents,  as  defined  in 
ARCSL-TR-77049^  was  modified  to  output  the  effective  time,  tg,  of  dosage  accumulation  as  a 
function  of  downwind  distance.  Data  were  generated  for  a variety  of  conditions  and  fitted  to  a 
general  equation. 


where 


^e 


0.005 


j^O.9294 


u 


(14) 


tg  is  the  effective  time  of  dosage  accumulation  (min) 

0.9294  is  the  slope  of  the  sigma-x  curve 

Given  this  effective  time  of  exposure,  one  can  then  consider  the  time  correction 
factor  for  the  dosage^  to  compute  the  effective  dosage  at  any  distance. 

M = 0.827  t>2min  (15) 

where 

M is  the  factor  by  which  the  required  dosage  is  increased 
t is  the  time  of  dosage  accumulation  (min) 


3.1.1  Two-Minute  Corrected  Dosage:  Segment  1,  x < X| 

By  setting  t = tg  and  multiplying  D by  the  factor  M (combining  equations  7,  14,  and 
1 5),  the  following  equation  can  be  used  to  calculate  distance  in  segment  1 . This  equation  is  also 
subject  to  the  condition  that 


13 


X > [400  u] 


1.076 


which  represents  the  distance  at  which  the  time  of  dosage  accumulation  will  exceed  2 minutes  as 
calculated  from  equation  14. 


X 


0 

36.51  ayia2iDu°-'^26 


1/a + /J  + 0.255 

for  Xj  > X > Xj 


(16) 


where  xj  = (400  u)  1-076 

For  X < Xj  and  xj,  use  equation  7 


3.1.2  Two-Minute  Corrected  Dosage:  Segment  3,  x > X2 

By  combining  equations  9,  14,  and  15,  the  2-minute  correction  in  segment  3 can  be 
computed  subject  to  the  same  restriction  in  tg. 


X 


0 

29.13  Uyi  Du0-'^26 


1/a + 0.255 

X > X2  and  xj 


(17) 


3.1.3  Two-Minute  Correction  of  Dosage:  Segment  2,  Xj  < x < X2 

Following  the  same  approach  outlined  above,  when  x > xj  , 


Dl  = 


36.51  CTyj  CT^i 


y0.726x  a + /J  + 0.255 


(18) 


D2  = 


29.13  Oyi  u 


0.726  ^ a + 0.255 


(19) 


The  distance  is  then  computed  using  equations  12  and  13.  If  x < x-p,  then 
equations  10  and  1 1 are  used  instead  of  equations  18  and  19. 

3.2  Semicontinuous  Source  with  2-Minute  Correction.  The  special  computer  program 
referred  to  in  section  3.1  was  executed  for  a range  of  release  times  to  provide  estimates  of  tg  for  the 
semicontinuous  release.  These  data  were  fitted  to  a generalized  model  as  follows: 

(20) 

L -I 

where  t^  is  the  source  release  time  (min). 
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. For  small  values  of  x,  equation  20  takes  on  a constant  value  based  on  the  value  of  tg. 

At  large  values  of  x,  the  values  are  the  same  as  those  obtained  from  equation  14.  In  the  transition 
between  these  two  extremes,  the  value  of  tg  is  underestimated  slightly  which  in  turn  provides  a 
safe-sided  estimate  of  hazard  distance. 


3.2.1  Semicontinuous  Source  with  2-Minute  Correction.  Segment  1 

Following  the  approach  outlined  in  section  3. 1 . 1 , the  equation  for  the  hazard 
distance  to  dosage,  D,  from  source,  Q,  is  as  follows: 


.2  ^ 0.000025  1.8588 

0.281  t^  + X 

^ 2 


0.137 


rOi  + P _ 


155.9  Oyj  u D 


xi>x>xjs  (21) 


where 


XTS 


O 0-,  0.538 
(4- 0.281  t2)u2 
5_ 


0.000025 


(22) 


For  values  of  x < x j and  Xj,  equation  7 is  used. 

It  is  noted  that  equation  21  is  not  solved  for  the  downwind  distance,  x.  Due  to  the 
form  of  the  equation,  it  was  decided  to  arrive  at  a value  of  x by  successive  approximation.  This 
approach  is  practical  on  a programmable  electronic  calculator  of  the  pocket-size  variety. 
(Experience  has  shown  that  some  solutions  can  take  a minute  or  more,  but  it’s  automatic  and  it 
does  get  there.) 


following: 


3.2.2  Semicontinuous  Source  with  2-Minute  Correction.  Segment  3 

Distances  that  fall  in  segment  3 are  computed  in  a similar  manner  using  the 
-1  0.137 


0.281 1^  + X 

^ ii2 


1.8588 


x“  = 


124.4  Oyj  Hj^  D u 


X > X2  and  xjg  (23) 


3.2.3  Semicontinuous  Source  with  2-Minute  Correction:  Segment  2,  x > xjg, 

The  values  of  Dj  and  D2  are  computed  in  segment  2 from  the  following  equations. 

^ r “1-0.137 

D.  = ^ 0.281  t2+M22225  ^ 1.8588  (24^ 

‘ 155.9ayia,iux“  + /3  [ ^ u2  * J 


D^ 


124.4  Oyj  Hj^  u x“ 


0.281  tf  + 


0.000025  1.8588 


] 


-0.137 


(25) 
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4. 


ACCURACY  OF  THE  SIMPLIFIED  METHODOLOGY 


The  methodology  has  been  tested  against  that  programmed  in  ARCSL-TR-77049 
over  a wide  range  of  input  parameters.  The  estimated  distances  are  always  found  to  equal  or  exceed 
the  reference  values  as  planned.  These  differences  are  generally  below  3%  with  some  few  ranging  as 
high- as  5%. 

5.  COMPUTATIONS  WITH  THE  SIMPLIFIED  METHODOLOGY 

With  the  exception  of  equations  21  and  23,  all  of  these  equations  can  readily  be 
solved  on  a pocket  calculator  which  has  a fractional  power  function  (y^).  The  only  complication  is 
the  determination  of  the  proper  segment  and,  thus,  the  selection  of  the  proper  equation  for  x.  The 
complete  logic  for  making  this  selection  is  shown  in  the  figure. 

The  procedure  starts  with  a total  dosage  solution  and  then  repeats  the  procedure 
using  the  2-minute-corrected  equations  when  needed.  The  figure  represents  a computer  approach 
which  is  complete  but  nonintelligent.  Once  the  user  has  a feel  for  the  distances  involved,  the  total 
dosage  solution  or  any  unneeded  segment  solution  may  be  by-passed  on  the  values  of  xj,  X2,  and 
X7  as  long  as  the  restrictions  on  the  basic  equations  are  complied  with. 

The  procedure  for  selecting  the  proper  segment  has  been  explained  in  the 
development  of  the  methodology.  The  Xj  restriction  is  somewhat  more  subtle.  The  time  of  dosage 
accumulation  must  exceed  2 minutes  before  the  2-minute  correction  is  applied.  Thus,  for  distances 
less  than  Xj  (or  xj§),  the  total  dosage  and  2-minute-corrected  solutions  are  the  same. 

6.  FORTRAN  PROGRAMS 

In  order  to  test  the  simplified  methodology  and  the  logic  in  the  figure,  two  Fortran 
computer  programs  were  written.  These  programs  contain  42  and  56  statements  as  compared  to  334 
for  subroutine  DOSDIS.^  One  provides  solutions  for  instantaneous  releases  and  the  other  for 
semicontinuous  releases.  Each  program  has  the  option  of  solving  either  the  total  dosage  or 
2-minute-corrected  methodology.  These  programs  are  very  similar  and  one  program  could  provide 
all  four  solutions  if  provisions  were  made  to  select  the  proper  reference  sigma  Y values  for 
instantaneous  and  continuous  solutions.  The  two  Fortran  programs  are  listed  in  appendix  A,  where 
the  inputs  are  defined. 

7.  POCKET  CALCULATOR  PROGRAMS 

In  the  first  attempts  to  code  this  methodology  on  a programmable  pocket 
calculator,  a two-segment*  approximation  to  the  mixing  layer  was  used.  A segment  of  the 
instantaneous  model  with  2-minute  correction  was  then  programmed  and  finally  the 
semicontinuous  with  2-minute  correction  using  an  iterative  approach  to  solve  equation  21  was 
tested.  Each  of  these  was  a program  of  about  a hundred  steps  and  was  written  for  the  Texas 
Instruments  Model  52. 


*The  intersection  of  the  Gaussian  and  Box  Models  separating  the  segments. 
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COMPUTE 
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UTE  X 

I 
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SET  X, 


> ANSWER 

T.  y 


Figure.  Logic  for  Determining  Proper  Segment 
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Once  the  complete  logic  was  tested  with  a Fortran  program,  this  approach  was- 
coded  for  the  Texas  Instruments  Model  59.  With  the  increased  memory  of  the  TI  59,  it  became 
possible  to  store  the  stability-dependent  parameters  (ayj,  a,  P,  Cj,  C2)  in  memory  and  select 
the  proper  set  by  indirect  addressing.  The  larger  program  size  made  it  possible  to  store  the 
three-segment  logic  as  defined  in  the  figure  including  the  2-minute  corrections.  Again,  two  programs 
are  used,  one  for  instantaneous  and  one  for  semicontinuous  sources.  These  programs  are  discussed 
in  appendix  B. 

8.  LIMITATIONS  OF  THE  SIMPLIFIED  METHODOLOGY 

The  simplified  methodology  treats  only  point  releases  at  ground  level  as  calculated 
from  a Gaussian  diffusion  model  with  vertical  boundaries.  The  2-minute  correction  as  defined  in 
ORG  40^  and  Technical  Paper  No.  10^  is  included  and  may  be  applied  to  G and  V dosages  as 
required.  The  simplified  methodology  does  not  include  height  of  release  or  source  cloud 
dimensions. 


Special  models  such  as  that  recommended  for  explosive  dissemination  of  HD  and  the 
combined  inhalation  and  disposition  of  VX  are  not  treated.  These  models  should  lend  themselves  to 
a similar  treatment.  However,  to  avoid  further  delay,  it  was  decided  to  complete  the  current  report 
at  this  time  and  follow  it  with  further  extensions  as  those  become  available. 

This  simplified  methodology  is  intended  as  a short-cut  approximation  and  the 
program  reported  in  ARCSL-TR-77049^  should  still  be  used  for  reference  and  for  special 
applications  not  treated  in  this  report. 
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APPENDIX  A 


FORTRAN  PROGRAMS 

Program  SINST  estimates  the  distance  to  a specified  dosage  of  interest  for  either 
total  dosage  or  2-minute-corrected  dosage. 

% 

When  program  SINST  is  executed,  it  will  request  the  following  information: 

INPUT:  Q,  u,  D,  HML,  IS,  IT 

1 

< where 

! Q is  source  strength  (mg) 

I 

i 

/ u is  windspeed  (m/sec) 

D is  dosage  (mg-min/m^) 

HML  is  height  of  the  mixing  layer  (m) 

IS  is  one  of  digits  1 through  6 corresponding  to  stability  categories  A through  F 
IT  is  an  indicator,  IT  = 0 for  total  dosage 

IT  = 1 for  2-minute  correction 

Variable  format  is  followed.  All  real  numbers  are  written  with  a decimal  point  and 
fields  are  separated  by  commas. 

The  program  outputs  Xj , X2,  and  x, 

where 

X j is  the  distance  at  which  reflection  from  the  inversion  cap  begins 
X2  is  the  distance  at  which  vertical  mixing  is  uniform 
X is  the  distance  at  which  the  dosage  of  interest  occurs. 
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WRK(1).SINST 

1 C SH0RT  T5?TAL  DHSAGE  (IT:0)/SH3RT  2-MIN  CRT  ( I T:  1) /S  A0/CGW 

2 DIMENSI0N  SYl ( 6) , ALF ( 6) , SZ 1 ( 6 ) , RTA (6 ) ,C 1 ( 6 ) ,C2 ( 6 ) 

3 DATA  SYl /. 09,  .0633,  .048,  .0634,  .0754,  ,07<56/ 

4 DATA  SZl /. 0222, . 1 1 ,. 1 19, .0898, .0879, .0791/ 

5 DATA  ALF/1 . , 1 . , 1 . , .9,  .8,  .7/ 

6 DATA  BTA/1 .4, 1 ., .9, .85, .8, .75/ 

7 DATA  C 1 /10.5 ,5 .5,6. 13,9.49, 1 1 .2, 15.3/ 

8 DATA  02/15.9,9.55,11.2,18.1,21.9,31.1/ 

9 3 PRINT  101 

10  101  FORMAT  ('INPUT:  Q, U , D, HML, IS , IT ' ) 

11  READ  100,Q,U,D,HML, IS,  IT 

12  APB:ALF(IS)+BTA(IS) 

13  HRB  = HML**(  1 ./BTA(IS) ) 

14  ALFI:ALF(IS) 

15  XT:1.E36 

16  IF  (IT  .EQ.  1)  XTr (400.*U)**1 .076 

17  XI :C1 (IS)*HRB 

18  X2=C2(IS)*HRB 

19  QYZDU=Q/(188.4*SY1 ( IS) *SZ 1 ( IS ) *D*U ) 

20  QYHDU=QYZDU*1 .253+SZl (IS) /HML 

21  5 X=QYZDU**( 1 ./APB) 

22  IF  (X  .LT.  XI ) GO  T3  4 

23  X=GYHDU**( 1 ./ALFI) 

24  IF  (X  .LT.  X2)  GO  TO  2 

25  4 IF  (X  .LT.  XT)  G0  TO  1 

26  U£=U** .274/. 1937 

27  QYZDU:QYZDU*UE 

28  GYHDU  = 0YHDU*UE 

29  APB:APB+.255 

30  ALFUALFI+.255 

31  XT=1.E36 

32  G3  TO  5 

33  2 D1 =QYZDU/X1**APB 

34  D2:QYHDU/X2**ALFI 

35  SrL3G(X2/Xl )/L3G(D2/Dl ) 

36  X=X1/D1**S 

37  GO  TO  4 

38  1 PRINT  200, XI ,X2,X 

39  GO  TO  3 

40  100  FORMAT  ( ) 

41  200  F0RMAT(23X, ' X 1 = ' , 1 PE  1 0 . 3 , ' X2=',E10.3,'  X=',E10.3,'  M') 


42 


END 


Appendix  A 
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Program  SSCS  estimates  the  distance  to  a specified  dosage  of  interest  from  a 
semicontinuous  (uniform  release  for  a finite  time)  source.  The  estimate  is  made  for  either  total 
dosage  or  2-minute  correction.  When  program  SSCS  is  executed,  the  following  input  is  requested: 

INPUT:  Q,  TS,  u,  D,  HML,  IS,  IT 


where 


TS  is  the  release  time  (min) 

Other  inputs  and  outputs  are  the  same  as  for  program  SINST. 


Appendix  A 
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WRK(l) .SSCS 
1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


C SH;1RT  SEMI-C0NT  (IT=0)  W/2-MIN  C0RR  ( I T=  I ) /SAO/CGW 
DIMENSION  SYl (6),ALF(6),SZ1 (6) ,BTA(6) ,CI (6) ,C2(6) 

DATA  S Y i / . 27 , . 1 899 , . 1 25 , . 1 268 , . 1 508 , . I 5 92/ 

DATA  SZ 1/. 0222,. n I 19,. 0898,. 0879,. 0 791/ 

DATA  ALF/1.,1.,1.,.9,.8,.7/ 

DATA  BTA/1.4,1.,.9,.85,.8,.75/ 

DATA  Cl /lO. 5, 5. 5,6. 13, 9. 49,1 1.2, 15.3/ 

DATA  02/15.9,9.55, 11.2,18.1,21.9,31.1/ 

PRINT  101 

101  F0RMAT(’  INPUT:  Q, TS, U,D,HML, IS, IT* ,6X, 'XI ' ,8X, ’X2’ ,8X, 


1 1 
12 
13 
1 4 

15 

16 
1 7 
18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
33 

39 

40 

41 

42 

43 

44 

45 

46 

47 
43 

49 

50 

51 

52 

53 

54 

55 

56 


> 


3 READ  100,Q,TS,U,D,HML,IS, IT 
APB:ALF(IS)+BTA(IS) 

HRB:HML**(1 ./BTA(IS)) 

ALFUALFdS) 

ET=.281+TS+TS 

RU:.000025/(U*U) 

X1=C1 (IS)*HRB 
X2=C2(IS)*HRB 

QYZDU=Q/( 188.4+SYl (IS)^SZl (IS)*D*U) 

QYHDU:QYZDU*1 .253+SZl (IS)/HML 
IF  (IT  .EQ.  0)  G0  T0  5 
XT  = 0. 

IF  (TS  .GT.  3.77)  G0  T0  6 
XT=((4.-ET)/RU)**.538 

5 X=QYZDU*+( 1 ./APB) 

IF  (X  .LT.  XI)  G0  T0  4 
X=QYHDU+*(1 ./ALFI) 

IF  (X  .LT.  X2)  G0  T0  2 

4 IF  (X  .LT.  XT  .0R.  IT  .EG.  0)  GO  TO  1 

6 CALL  XE(QYZDU,APB) 

IF  (X  .LT.  XI)  G0  T3  1 
CALL  XE(Q  YHDU,ALFI) 

IF  (X  .GT.  X2)  G0  T3  1 

QYZDU:QYZDU/(.827*(ET+RU*X1** 1 .8533)**. 137) 
QYHDU:QYHDU/(.827*(ET+RU*X2**1 .8533 )**. 137) 

XT=1 .E36 

2 D1=QYZDU/X1**APB 
D2:QYHDU/X2**ALFI 
S=L0G(X2/X1 )/L0G(D2/Dl ) 

X=X1/D1**S 
G0  T0  4 

1 PRINT  200, XI  ,X2,X 
G0  T0  3 
100  FORMAT  ( ) 

200  FORMAT  (25X,3F10.0) 

SUBROUTINE  XE(FQ,P) 

E=FQ/.827 

X=0. 

DX=1 .E5 

7 X=X+DX 

IF  (((ET+RU*X**1 .8583)**. 137*X**P)  .LT.  E)  GO  T0  7 

X=X-DX 

DX=DX/10. 

IF  (DX  .GE.  1 .)  G0  T0  7 

RETURN 

END 


X ( M)  * 


Appendix  A 
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APPENDIX  B 


POCKET  CALCULATOR  PROGRAMS 


The  programs  listed  in  this  appendix  were  written  for  the  Texas  Instruments 
pocket  calculator  Model  59.  The  programs  could  readily  be  adapted  to  any  calculator  of  similar 
capability.  These  programs  use  36  words  of  memory  to  store  the  stability-dependent  parameters 
and  6 words  to  store  the  indirect  address  references.  This  approach  has  the  advantage  that  the 
stability  is  simply  selected  by  the  user  and  the  parameters  do  not  have  to  be  keyed  in  with  each 
execution.  For  program  SINST,  1 2 additional  words  are  used  for  intermediate  storage,  and 
program  SSCS  uses  1 8 additional  words.  Program  SINST  occupies  260  instruction  cells,  and  SSCS 
occupies  419  instruction  cells. 

The  pocket  calculator  programs  follow  the  logic  diagram  (figure  in  text)  and  the 
Fortran  programs  listed  in  appendix  A very  closely.  The  user  should  note  that  the  same  input  is 
used  but  the  order  is  different. 

For  the  TI  59,  either  program  is  loaded  from  two  magnetic  cards.  (Systems 
Assessment  Office  will  write  three  programs  on  your  cards  if  requested.)  The  program  is  then 
initiated  by  pressing  label  A.  The  following  tables  wUl  indicate  the  order  of  entry  for  each 
program. 


TABLE  B-1.  DATA  ENTRY  PROGRAM  SINST 

Label  A IS  is  one  of  digits  1 through  6 to  indicate  stability,  A through  F 

Label  B HML  is  height  of  the  mixing  layer  (m) 

Label  C Q is  quantity  airborne  (mg) 
u is  windspeed  (m/sec) 

D is  dosage  of  interest  (mg-min/m^) 

IT  =0  for  total  dosage 

= 1 for  2-minute  correction 
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TABLE  B-2.  DATA  ENTRY  PROGRAM  SSCS 


Label  A 

IS  is  one  of  digits  I through  6 to  indicate  stability,  A through  F 

Label  B 

HML  is  height  of  the  mixing  layer  (m) 
tg  is  release  time  (min) 

Label  C 

Q is  quantity  airborne  (mg) 
u is  windspeed  (m/sec) 

D is  dosage  of  interest  (mg-min/m^) 
IT  =0  for  total  dosage 

= 1 for  2*minute  correction 

The  program  will  progress  from  Label  A to  B to  C automatically  and  return  to 
Label  C with  the  estimated  distance.  If  the  user  does  not  want  to  restart  with  Q,  he  may  restart 
at  Label  B or  A. 
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DISTRIBUTION  LIST  4 


Names 

CHEMICAL  SYSTEMS  LABORATORY 

SAFETY  OFFICE 
Attn:  DRDAR-CLF 

PLANS  & PROGRAMS  OFFICE 
Attn:  DRDAR-CLR-L 
Attn:  DRDAR-CLR-T 

AUTHOR’S  COPIES:  Systems  Assessments  Office 

BIOMEDICAL  LABORATORY 
Attn:  DRDAR-CLL-B 
Attn:  DRDAR-CLL-M 
Attn:  DRDAR-CLL-V 

CB  DETECTION  & ALARMS  DIVISION 
Attn:  DRDAR-CLC-C 
Attn:  DRDAR-CLC-E 

DEVELOPMENTAL  SUPPORT  DIVISION 
Attn:  DRDAR-CU-R 
Attn:  DRDAR-CU-L 
Attn:  DRDAR-CU-M 
Attn:  DRDAR-CU-P 

ENVIRONMENTAL  TECHNOLOGY  DIVISION 
Attn:  DRDAR-CLT-D 

MUNITIONS  DIVISION 
Attn:  DRDAR-CLN 


PHYSICAL  PROTECTION  DIVISION 


Attn: 

DRDAR-CLW 

Attn: 

DRDAR-CLW-C 

Attn: 

DRDAR-CLW-P 

Attn: 

DRDAR-CLW-E 

RESEARCH  DIVISION 

Attn: 

DRDAR-CLB 

Attn: 

DRDAR-CLB-B 

Attn: 

DRDAR-CLBC 

Attn: 

DRDAR-CLB-P 

Attn: 

DRDAR-CLB-R 

Attn: 

DRDAR-CLB-T 

Attn: 

DRDAR-CLB-TE 

SYSTEMS  ASSESSMENTS  OFFICE 

Attn: 

DRDAR-CLY-A 

Attn: 

DRDAR-CLY-R 

DEPARTMENT  OF  DEFENSE 

Administrator 

Defense  Documentation  Center 

Attn:  Accessions  Division  (DDC-TC) 
Cameron  Station 
Alexandria,  VA  22314 


Names  Copies 

Director 

Defense  Intelligence  Agency 

Attn:  DB4G1  1 

Washington,  DC  20301 

DEPARTMENT  OF  THE  ARMY 

HQDA  (DAMO-SSQ  I 

WASH  DC  20310 

Director 

Defense  Qvil  Preparedness  Agency 

Attn:  PO(DC)  1 

Washington,  DC  20301 

CINCUSAREUR 

Attn:  AEAGC-RSI  1 

APO  New  York  09403 

Deputy  Chief  of  Staff  for  Research, 

Development  & Acquisition 

Attn:  DAMA-CSM-CM  1 

Attn:  DAMA-ARZ-D  1 

Washington,  DC  20310 

US  Army  Research  and  Standardization 

Group  (Europe)  1 

Attn:  Chief,  Chemistry  Branch 
Box  65,  FPO  New  York  09510 

HQDA  (DAMI-FIT)  1 

WASH,  DC  20310 

Commander 

HQ  US  Army  Medical  Command,  Europe 

Attn:  AEMPM  1 

APO  New  York  09403 


US  ARMY  HEALTH  SERVICE  COMMAND 

Superintendent 
Academy  of  Health  Sciences 
US  Army 

Attn:  HSA-CDC  1 

Attn:  HSA-IHE  1 

Fort  Sam  Houston,  TX  78234 

US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 

US  Army  Materiel  Development  and  Readiness  Command 
Attn:  DRCLDC  1 

Attn:  DRCSF-P  1 

5001  Eisenhower  Ave 
Alexandria,  VA  22333 


Copies 

I 

4 

1 

40 

1 

1 

1 
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1 

3 
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1 

I 

4 

2 

2 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

4 

1 

2 


35 


DISTRIBUTION  LIST  4 (Contd) 


Names  Copies 

Office  of  the  Project  Manager  for  Chemical  Demilitarization 
and  Installation  Restoration 

Attn:  DRCPM-DR-T  2 

Aberdeen  Proving  Ground,  MD  21010 

Human  Engineering  Laboratory  HFE  Detachment 

Attn:  DRXHE-EA  1 

Building  E3220 
APG'Edgewood  Area 

Commander 

US  Army  Foreign  Science  & Technology  Center 

Attn:  DRXST-CX2  1 

220  Seventh  St.,  NE 
Charlottesville,  VA  22901 

Commander 

US  Army  Missile  Research  and  Development  Command 
Redstone  Scientific  Information  Center 

Attn:  DRDMI-TBD  1 

Redstone  Arsenal,  AL  35809 

Director 

DARCOM  Field  Safety  Activity 

Attn:  DRXOS-C  1 

Charlestown,  IN  47111 

Commander 

US  Army  Materiel  Development  and  Readiness  Command 
Installations  and  Services  Activity 

Attn:  DRCIS-RMU  1 

Rock  Island,  IL  61202 

US  ARMY  ARMAMENT  RESEARCH  AND 
DEVELOPMENT  COMMAND 

Commander 

US  Army  Armament  Research  and  Development  Command 


Attn:  DRDAR-AC  1 

Attn:  DRDAR-LCEtM  1 

Attn:  DRDAR-LCH  1 

Attn:  DRDAR-LCN  1 

Attn:  DRDAR-LCU  1 

Attn:  DRDAR-SER  1 

Attn:  DRDAR-TSS  2 

Attn:  DRCPM-CAWS-WP  1 

Dover,  NJ  07801 

Director 

Ballistic  Research  Laboratory 

Attn:  DRDAR-TSD-S  1 

Building  328 
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